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The first measurement of D*^ meson diffractive pfioto- 
production cross sections fias been performed witfi tfie ZEUS 
detector at the HERA ep collider, using an integrated luminosity 
of 38pb~^. The measurement has been performed for photon- 
proton center-of-mass energies in the range 130 < W < 280 GeV 
and photon virtualities < 1 GeV^. D*^ mesons have been 
reconstructed with > 2 GeV and — 1.5 < r]^* < 1.5 from the 
decay channel L>*+ D^7r+ with K-7T+ (+c.c.). The 

diffractive component has been selected with 0.001 < Xp < 0.018. 
The measured cross section in this kinematic range is: 
crfJSe'D^Xp' = 0-74 ± 0.21(stot)^!!:?^(s7/st.) ± {).lQ{p.diss.) nb (ZEUS 
preliminary). Measured integrated and differential cross sections have 
been compared to theoretical expectations. 
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1 Introduction 



Charm production processes are the ones which proceed mainly through 
gluon-initiated hard subprocesses and are perturbatively calculable. Thus, 
the diffractive production of charmed mesons can provide new tests of the 
partonic structure of diffractive interactions, in particular of their gluon 
component. 

During the years of the HERA collider operation, integrated and dif- 
ferential cross sections for inclusive charm production were measured in 
kinematic ranges where an effective signal separation from suppressed back- 
grounds could be achieved The measured cross sections for photopro- 
duction (PhP) and deep inelastic scattering (DIS) processes were compared 
with different next-to-leading (NLO) pQCD calculations. DIS data were 
found to be in good agreement with the calculations. The calculated PhP 
cross sections are lower than the measured ones, especially in the forward 
(proton) direction. 

As for the diffractive charm production, there are only preliminary 
results on diffractive dissociation in DIS, measured with D*^ mesons 
Here we present preliminary results on measurements by the ZEUS Col- 
laboration of cross sections for diffractive photoproduction of D*^{2010) 
mesons]^ in the range of Pomeron fractional momentum 0.001 < Xp < 0.018 
at energies 130 < W < 280 GeV in the photon-proton center-of-mass frame 
and photon virtualities < 1 GeV^. D* mesons were reconstructed through 
the decay channel D*+ D^7r+ (K^7r+)7r+ (and c.c.) in the restricted 
kinematic region: p^* > 2 GeV and \r]^*\ < 1.5. Here p^* is the D* me- 
son transverse momentum and r/^* = -ln(tan(^/2)) is its pseudorapidity, 
defined in terms of the D* polar angle 6 with respect to the proton beam 
direction. 

The measurements were performed at the HERA collider with the 
ZEUS detector, a detailed description of which can be found elsewhere 
[jSj. The data were taken during 1996 and 1997, when HERA colhded 
positron and proton beams with energies of 27.5 GeV and 820 GeV, re- 
spectively. An integrated luminosity of 38 pb^^ was used for this measure- 

^ In the following, 13*^(2010) will be referred to simply as D* . 
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ment. Charged particles were measured in the central tracking detector 
(CTD) [0. To detect the scattered electron and to measure global energy 
values the uranium-scintillator sampling calorimeter (CAL) [[7| was used. 
The luminosity was determined from the rate of the bremsstrahlung pro- 
cess e~^p e^7j9, where the photon was measured by a lead scintillator 
calorimeter 



2 Kinematics of difFractive photoproduction 

We consider diffractive photoproduction in ep scattering at HERA: 

e(e) + p{p) e'(e') + X + p'{p'), 

where the four-momenta of particles are shown in brackets. The collision 
occurs at the squared positron-proton center-of-mass energy s = (e + p)'^, 
and photon virtuality = — g^, where q = e — e'. The squared photon- 
proton center-of-mass energy = {p + is defined for this reaction. 
One may consider that the interaction proceeds through a photon-Pomeron 
(iP) scattering: 

7(g) + iP(P.) ^ X, 

where Pp = p — p'. This process is described by the invariant mass Mx of 
the hadronic system X, produced by photon dissociation, and the fraction 
of the proton momentum 

Pjp-q Ml 
"~ p-q - 

carried away by the Pomeron. 

The variables VF, Mx and Xp were reconstructed from the final hadronic 
system, measured by energy flow objects (EFO) 0], made from tracks 
detected by the CTD and from energy deposits in the CAL cells. The 
Jacquet-Blondel formula W,^ = ^j2EpY,^{E - P,)^ [[TOj was used to re- 



construct W . Here Ep is the proton beam energy. The invariant mass of 
the diffractively produced system Mx was calculated with the formula 

Ml = (E. E,f - (E. P.f - (E. Pyf - (E. Pzf- Sums in both equa- 
tions run over energies Ei and momenta Pi of all EFOs. To select the 
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W range, VFjb was calculated with calorimeter cell deposits only so as to 
be consistent with the inclusive charm photoproduction analysis ^ . Mea- 
sured values were corrected to the true ones by factors, determined from 
MC simulations of diffraction as average ratios of reconstructed to gener- 
ated values. All variables were reconstructed to an accuracy of better than 
15%. 



3 Event selection and D* reconstruction 



Event selection and reconstruction procedures are described in 
details elsewhere |T|. Here a short description is given. 

Photoproduction events were selected by requiring that no scattered 
positron was identified in the CAL [jTT|] and the photon-proton center-of- 
mass energy W is between 130 and 280 GeV. Under these conditions, 
the photon virtuality is limited to values less than 1 GeV^. The cor- 
responding median was estimated from a Monte Carlo (MC) simu- 
lation to be about 3 x 10^^ GeV^. The mesons were reconstructed 
through the decay channel D* (D^ Ktt)tTs by combining candidates 
from charged tracks measured by the CTD. For the reconstruction, "right 
charge" track combinations, defined for [Kt() with two tracks of oppo- 
site charges and with a tt^ having the charge opposite to that of the 
K meson in the (i^7r), were accepted as long as the combination of in- 
variant masses AM = M{KTnis) — M{Kti) and M{Ktt) are within wide 
mass-windows around the nominal values of AM = M{D*) — M{D^) and 
M{D^) [|T2|. To determine the number of D* mesons in the signal, combi- 



natorial background was modelled by "wrong charge" track combinations 
and subtracted after normalization to the "right charge" distribution in 
the range 0.15 < AM < 0.17 GeV. "Wrong charge" combinations were 
defined for [Kn) with two tracks of the same charge and with a tt^ of the 
opposite charge. The measurements were performed in the pseudorapidity 
range — 1.5 < t]^* < 1.5, where the CTD acceptance is high. The kinematic 
region in was limited to 2 < < 8 GeV. 

The MC event samples used for this analysis were prepared with 
the RAPGAP 0, PYTHIA and HERWIG generators. Diffrac- 
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tive interactions were modeled in the framework of the resolved Pomeron 
model []TB| with /3{1 — (3) or the HI FIT2 []T7|] parametrisations for the 



initial partonic distributions in the Pomeron. Here (3 is the fraction of the 
Pomeron momentum carried by a parton, that couples to the Pomeron and 
participates in the hard interaction. The MRSG [|TB| and GRV-G HO [jT^ 



parametrisations were used for the proton and photon structure functions, 
respectively, when modelling non-diffractive interactions. The fragmenta- 
tion of the generated partons (parton shower evolution and hadronisation) 
was simulated according to the LUND model when using the RAPGAP 



or PYTHIA simulations. The HERWIG generator models the hadronisa- 
tion process with a cluster hadronisation model. The MC events were 
processed through the standard ZEUS detector and trigger simulation pro- 
grams and through the same event reconstruction package as was used for 
data processing. The shapes of MC and data distributions were found to 
be in reasonable agreement within statistical errors. 

Diffractive events were identified by a large rapidity gap (LRG) be- 
tween the scattered proton, which escaped detection through the beam 
pipe, and the hadronic system X, produced by the dissociated photon. 
The LRG events were searched for using the ijmax method, for which ijmax 
was defined as the pseudorapidity of the most forward EFO with energy 
greater than 400 MeV. Fig. 1 presents the r]„^ax distribution for all photo- 
produced D* mesons, reconstructed within the signal range 
0.143 < M{KTnTs) - M{Ktt) < 0.148 GeV and 1.80 < M{K7t) < 1.92 GeV 
after the combinatorial background subtraction. This distribution shows 
two structures. The plateau-like structure at ijmax ^ 2 is populated pre- 
dominantly by the LRG events, while the wide peak-like structure around 
Vmax ~ 3.5 originates from the non-diffractive events and has an exponen- 
tial fall-off towards lower values of ijmax- From a comparison between the 
data points and a sum of simulated diffractive and non-diffractive event 
distributions, normalised to the data, a cut-off of ijmax = 1-75 was chosen 
as a compromise between the magnitudes of the diffractive signal and the 
non-diffractive background. The non-diffractive background fractions for 
subtraction were estimated from the MC-to-data distribution ratios, using 
non-diffractive MC simulations. 

When using the ijmax method for diffractive event selection, one needs 
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to take into account the following properties of the method. The measure- 
ment of rapidities by the CAL is limited to the edge of the forward beam 
hole of the CAL. Thus the proton dissociative events, ep e'XN, can sat- 
isfy the requirement ijmax < 1-75 if the proton dissociative hadronic system 
N has invariant mass small enough to pass undetected through the forward 
beam pipe. It was found earlier that the proton dissociation contribution 
comprises 0.31 ± 0.15 Measured cross sections were corrected for 



this value. A cut in ijmax correlates with a range of accessible Xp values. 
ijmax < 1-75 restricts Xp < 0.018. In addition, limited acceptance restricts 
Xp > 0.001. 

After the above selection and the "wrong charge" background subtrac- 
tion a signal of 56 ± 10 diffractively photoproduced D* mesons was found 
in the AM distribution (Fig. 2). 



4 Cross sections 



The inclusive D* production cross section is given by: 

ATcorr 
D* 



0'ep^D*X — -p, 5" , 

where N'^I'^ is the number of observed D* mesons corrected for the accep- 
tance, jC = 38.0 ± 0.6 pb~^ is the integrated luminosity and 5£,*^po^^^)^ = 



0.0263 ± 0.0010 is the combined D* {D^ K+tt-) tt^ decay branching 
ratio []T^. Acceptance corrections were calculated using the RAPGAP MC 
sample. 

The total D* diffractive photoproduction cross section in the kine- 
matic region < i QeV^, 130 < W < 280 GeV, p^* > 2 GeV, \r]^*\ < 1.5 
and 0.001 < Xp < 0.018 was measured to be 

%^e'i?*xp' = 0-74±0.21(stot)^g:?^(s7/st)±0.16(j9.dzss.) nb (ZEUS pre- 
liminary). The last error is due to the uncertainty in the proton dissociative 
background subtraction. Other sources of systematic uncertainties due to 
analysis and detector features were studied and their effect on the cross 
section was estimated. The largest contributions to the systematic error 
came from the CAL energy scale uncertainty (^4^8^^%), the signal determi- 
nation procedure (^145%), the selection of diffractive events (^8o^%) and 
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the acceptance correction calculations (^^5 9%). The overall normalisation 
uncertainties due to the error in the luminosity value (±1.7%) and in the 
D* and decay branchings (±3.8%) were not included in the system- 
atic error quoted above. All of the systematic uncertainties were added in 
quadrature to determine the overall systematic uncertainty of ^^^\%. The 
summation of the systematic uncertainties was also performed for each bin 
of the differential distributions. 

The measured D* diffractive photoproduction cross section, while only 
a fraction of the total diffractive contribution, amounts ~ 4% of the inclu- 
sive D* photoproduction cross section, (Jep^D*x = 18.9±1.2(stat)]^Q8(s?/st) 
nb measured in the same kinematic range. This fraction indicates that 
diffractive charm production is not suppressed as much as some early mod- 
els predicted . 



Measurements were compared to resolved Pomeron model expecta- 
tions [M , calculated with the RAPGAP Monte Carlo program in the same 



kinematic region. Partonic distributions in the Pomeron were parametrised 
by the fit to the HERA data [|T7|, performed by HI Collaboration (HI 
FIT2). Only the BGF mechanism of charm production was accounted for. 
The leading order RAPGAP Monte Carlo, with the HI FIT2 Pomeron 
parametrisation, predicts 1.42 nb for the diffractive photoproduction 
cross section in the same kinematic range [|ZB|] . 

Differential cross sections for da /dp^\ da/di]^* , da/dMx and da/dxp 
are presented in Figs. 3-6. All of the above mentioned systematic uncer- 
tainies were added in quadratures with statistical errors (inner error bars) 
in each bin to calculate the total error (outer error bars), both of which 
are shown in Figs. 3-6. 

The measured differential cross sections (Figs. 3-6), when compared 
to the ones of the resolved Pomeron model calculated with the RAPGAP 
MC program show reasonable agreement in shape with the theoretical ex- 
pectations considering the measurement errors, da/dp^* agrees well and 
the other three distributions are shifted somewhat to larger values with 
respect to the predictions. 
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5 Summary and conclusions 

The first measurement of diffractive D* photoproduction has been 
performed with the ZEUS detector at HERA with a luminosity of 38pb~^. 
The preliminary results are reported here. The total D* diffractive photo- 
production cross section in the kinematic region (^^ < 1 GeV^, 130 <W < 280 
GeV, py* > 2 GeV, \r]^*\ < 1.5 and 0.001 <X]p < 0.018 is measured to be 
(TfJSe'D*Xp' = 0.74±0.21(stot.)lg:g(syst.)±0.16(p.t/iss.) nb (ZEUS pre- 
liminary) . The leading order calculations in the framework of the resolved 
Pomeron model predict 1.42 nb for this cross section. The differential 
cross section shapes for da/dp^*, da/drj^*, dcr/dMx and da/dxp show 
reasonable agreement with the resolved Pomeron model considering the 
measurement errors. 
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Figure 1: Comparison of the measured rj^ax distribution (dots) with the sum 
of the diffractive and non-difjractive MC distributions (histograms) for events with 
D* mesons. D* candidates were selected in the kinematic region < 1 CeV"^ , 
130 < < 280 GeV, p^' > 2 CeV and \r]^*\ < 1.5. Sum distribution of the 
diffractive resolved Pomeron RAP CAP MC (dotted histogram) and non- diffractive 
MC (dashed histogram) events were normalised to have the same area as the data 
distribution. 
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Figure 2: The AM distribution for the D* diffractive photoproduction reac- 
tion with D* KTi)ns for < i GeV^ , 130 < < 280 GeV and 
0.001 < Xp < 0.018. The kinematic range of measurements is p!^* > 2 GeV and 
|?7^*| < 1.5. The dots are for the right charge combinations, and the dashed his- 
togram is for the wrong charge combinations from the signal region ( 1.80 - 1.92 
GeV ). The full line is the result of a fit to a sum of a Gaussian and the functional 
form A(AM - m^)^ . 
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Figure 3: Differential cross section da/dpJ^* (dots) for the diffractive pho- 
toproduction reaction ep e'D*Xp' for < I GeV^, 130 <W< 280 GeV and 
0.001 < Xjp < 0.18. The kinematic range of measurements is pf^* > 2 GeV and 
\ri^*\ < 1.5. The inner bars show the statistical errors, and the outer bars cor- 
respond to the statistical and systematic errors, added in quadrature. The data 
are compared with the distributions of the RAPGAP MG diffractive events, simu- 
lated in the framework of the resolved Pomeron model with the HI FIT2 Pomeron 
parametrization ( histogram). The MG distribution has been normalized to have 
the same area as the data distribution. 
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Figure 4: Differential cross section da/drj^* (dots) for the diffractive pho- 
toproduction reaction ep e'D*Xp' for < I GeV"^, 130 <W< 280 GeV and 
0.001 < Xjp < 0.18. The kinematic range of measurements is pf^* > 2 GeV and 
\ri^*\ < 1.5. The inner bars show the statistical errors, and the outer bars cor- 
respond to thestatistical and systematic errors, added in quadrature. The data 
are compared with the distributions of the RAPGAP MG diffractive events, simu- 
lated in the framework of the resolved Pomeron model with the HI FIT2 Pomeron 
parametrization ( histogram). The MG distribution has been normalized to have 
the same area as the data distribution. 
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Figure 5: Differential cross section da/dMx (dots) for the diffractive pho- 
toproduction reaction ep — > e'D*Xp' for < 1 GeV'^, 130 <W< 280 GeV and 
0.001 < Xjp < 0.18. The kinematic range of measurements is pf^* > 2 GeV and 
\ri^*\ < 1.5. The inner bars show the statistical errors, and the outer bars cor- 
respond to thestatistical and systematic errors, added in quadrature. The data 
are compared with the distributions of the RAPGAP MG diffractive events, simu- 
lated in the framework of the resolved Pomeron model with the HI FIT2 Pomeron 
parametrization ( histogram). The MG distribution has been normalized to have 
the same area as the data distribution. 
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Figure 6: Differential cross section da/dxp (dots) for the diffractive pho- 
toproduction reaction ep e'D*Xp' for < I GeV"^, 130 <W< 280 GeV and 



0.001 < Xp < 0.18. The kinematic range of measurements is pij^ > 2 GeV and 
\ri^*\ < 1.5. The inner bars show the statistical errors, and the outer bars cor- 
respond to the statistical and systematic errors, added in quadrature. The data 
are compared with the distributions of the RAPGAP MG diffractive events, simu- 
lated in the framework of the resolved Pomeron model with the HI FIT2 Pomeron 
parametrization ( histogram). The MG distribution has been normalized to have 
the same area as the data distribution. 
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